H> v * , . 
*Cl« r ' 






NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




REAL TIME SIMULATION AND CONTROL 
3000 TON SURFACE EFFECT SHIP 
WITH NEGATIVE DRAG CHARACTERISTICS 
IN SEA STATE 

by 

Lee Lewis Oliphant 
December 1980 

Thesis Advisor: A. Gerba 

Approved for public release; distribution unlimited 



T 197809 



security classification of this p > r,c o«*« 






REPORT DOCUMENTATION PAGE 


read instructions 

BEFORE COMPLETING FORM 


MPORT mumSC* 


2. GOVT ACCESSION MO. 


». RIClPUNT'l CATALOG NUMBER 


4 TITLE (And Subtttlo) 

Real Time Simulation and Control 3000 Ton 
Surface Effect Ship With Negative Drag 
Characteristics In Sea State 


5. TYPE OF REPORT A PCRtOO COVCREO 

Master's Thesis; 
December 1980 


i. PCRPORMI n g org. REPORT NUMIER 


7. AU T NO Afo) 

Lee Lewis Oliphant 


• . CONTRACT OR grant nomACAmj 


• PERFOPMING ORGANIZATION NAME ANO AODAdl 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM ELEMENT. PROJECT, task 
ARIA * WORK UNIT NUMICRS 


II CONTROLLING a»»ICl NAME ANO AOORE1J 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

December 1980 


11. NUMIEl of PAGES 

117 


1 >4 MONITORING AGENCY name k •OO^CS %ni dittoront from Controlling Oillco) 


is. security class, (ot (*<• roponj 

UNCLASSIFIED 


ls«. OECL ASSl F| C A TlON/ OOWNGRAOinG 

schedule 


I*. OISTRISuTlON STATEMENT (a 1 mia K apart) 

Approved for public release; distribution unlimited 


IT. Ol ST ftl EU TlON ST ATCMCN T fol !*• oAmtrmct on food J#t Block 20. II dit lor on 1 from Koport) 


IS SUPPLEMENTARY notes 


It. KEY WON OS ( Contlnv on rowtoo m*dm If nmcooooty And identity kf klock numkor) 

Surface Effect Ship; Real Time Simulation; Negative Drag 
Characteristic; Speed Control; Sea State Effect 


20 AISTRACT fContimto on rowoomm aldo 11 nocoooorr And Idontlfy by klock mmtbot) 

The model of a Surface Effect Ship was refined to include simpli- 
fied propulsion dynamics, negative drag characteristics, sea state 
effects and an autopilot for speed control. These design modifi- 
cations were introduced into a real time, man controlled simulation 
of a 3000 ton Surface Effect Ship (3K-SES) in 5 degrees of freedom 
(RTS5D) and results were compared with a Data Base Program (DBSIM5D) 
based on towing tank data scaled up to model the ship. 



DO , 1473 coition or 1 MOV «* IS OBSOLETE 

(Pa^e 1) S/N 0102*014 AA0 1 J SECURITY CLASSIFICATION OF THIS * AOl Ootm Bntorod) 



^cu^tv cl Ai>ir<c * t<on or n»n M*ffd 



Hardware and Software design changes were incorporated into 
the RTS5D model to provide a more accurate approximation of 
real time, a faster computer iteration time, and a broach 
condition warning if the operator exceeded certain thrust 
vectoring limits. 



DD Form 1473 




2 



**cu«t*v cuAMiriCATio* 0* thi* 



Approved for public release ; distribution unlimited 



Real Time Simulation and Control 
3000 Ton Surface Effect Ship 
With Negative Drag Characteristics 
In Sea State 



by 

. Lee Lewis Oliphant 
Lieutenant, Unites States Navy 
B.S., University of Texas, 1976 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 



from the 



NAVAL POSTGRADUATE SCHOOL 
December 1980 



ABSTRACT 



The model of a Surface Effect Ship was refined to include 
simplified propulsion dynamics, negative drag characteris- 
tics, sea state effects and an autopilot for speed control. 
These design modifications were introduced into a real time, 
man controlled simulation of a 3000 ton Surface Effect Ship 
(3K-SES) in 5 degrees of freedom (RTS5D) and results were 
compared with a Data Base Program (DBSIM5D) based on towing 
tank data scaled up to model the ship. 
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into the RTS5D model to provide a more accurate approxima- 
tion of real time, a faster computer iteration time, and a 
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NOMENCLATURE 



A ws Average sidewall wetted area, 
starboard side 



A Average sidewall wetted area, • 

wp port side 

Aj^ Added mass coefficient in roll 

force equation 

A^j Added mass coefficient in pitch 

force equation 



A - Average wetted sidewall area 
of the bow 



A wl 
A 2 2 



s 



c 

c 

c 

c 

F 

F 

F 

F 

F 

F 



DX (K) 

DY 

DZP 

DP 

sw 

ss 

sp 

1 

2 

3 



cr 

o 



Average wetted sidewall area 
of the stern 

Added mass coefficient in yaw 
force equation 

Sideslip angle 

Coefficients drag in x-direction 

Coefficient drag in y-direction 

Sidewall roll moment lumped 
parameter coefficient 

Bow pitch force lumped parameter 
coefficient 

Sidewall roll force 

Sidewall starboard buoyancy force 

Sidewall port buoyancy force 

Stern buoyancy force 

Bow seal pitching force 

Bow buoyancy force 

Gravitational acceleration 



ft 

ft 



2 

2 



ft-slug 



ft-slug 



ft 



2 



ft 



2 



slug s 2 lb J 
ft 1 

rad 

lb,s 2 lb /ft 2 

lb,s 2 lb /ft 2 
r m 

lb f s 2 /ft 2 



non- 

dimensional 

lbs 

lbs 

lbs 

lbs 

lbs 

lbs 



ft/s 



2 
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1^. Moment of inertia about x-axis 

Iy Moment of intertia about y-axis 

I- Moment of inertia about z-axis 

K Summation of moments about x-axis 

& sw Length of sidewall 

Actual draft of port sidewall 

I x-direction displacement of hull 

drag centroid 

Jlj Average draft of SES sidewall 

a l 

Average draft of bow seal 

l Length from center of gravity to 

x stern 

i Average draft of widewall 

X 1 

Pitch moment lever arm for bow 
sidewall buoyance force 

i- Pitch moment lever arm for bow 

seal force 

m Mass of the rigid ship 

M Summation of moments about y-axis 

N Summation of moments about z-axis 

P^ Plenum pressure 

p Lumped drag centroid point 

p' Lumped drag centroid point 

p Roll acceleration 



slug ft“ 

slug ft^ 

slug ft" 

lbs - ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

lbs - ft 
lbs - ft 
lbs - f t" 
non- 

dimensional 

non- 
dimens ional 

rad/ sec^ 
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q 


Pitch rate 




rad/sec 


q 


Pitch acceleration 




rad/sec 


r 


Yaw rate 




rad/sec 


r 


Yaw acceleration 




rad/ sec 


S 1 


Turning moment lever arm of no. 
engine 


1 


ft 


S 2 


Turning moment lever arm of no. 
engine 


2 


ft 


S 3 


Turning moment lever arm of no. 
engine 


3 


ft 


S 4 


Turning moment lever arm of no. 
engine 


4 


ft 


T 

7 


Total thrust magnitude on no . 1 
engine 




lbs 


H 

oo 


Total thrust magnitude on no. 2 
engine 




lbs 


T 

9 


Total thrust magnitude on no . 3 
engine 




lbs 


T 

A io 


Total thrust magnitude on no. 4 
engine 




lbs 


T 

forw 


Total forward thrust vector of 
effectors 




lbs 


T 

side 


Total side thrust vector of 
effectors 




lbs 


T 

yaw 


Total turning moment generated 
effectors 


by 


lbs ft 


u 


Velocity in x-direction (surge) 




ft/ sec 


u 


Acceleration in x-direction 




ft/ sec" 


V 


Velocity in y-direction (sway) 




ft/sec 


V 


Acceleration in y-direction 




ft/sec 


V 


Total velocity 




ft/ sec 
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w 

e 


Width of bow seal 




ft 


X 


Summation of forces in x 


-direction 


lbs 


X o 


X...,. coordinate of SES 
NAV 




ft 


X 

0 


X NAV velocity of SES 




ft/sec 


Y 


Summation of forces in y 


-direction 


lbs 


Y 

0 


Ynav coordinates of SES 




ft 


Y 

0 


Ynav velocity of SES 




f t/sec 


6 


Effector angle commanded 




rad 


5 7 


Effector angle of no. 1 


nozzle 


rad 


*8 


Effector angle of no. 2 


nozzle 


rad 


6 9 


Effector angle of no. 3 


nozzle 


rad 


5 10 


Effector angle of no. 4 


nozzle 


rad 


<P 


Heading angle of SES 




rad 




Roll angle of SES 




rad 


9 


Pitch angle of SES 




rad 


P 


Density of sea water 




lbm/f t 
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I. INTRODUCTION 



Rising fuel costs, critical manpower shortages and an 
increasing need for rapid deployment of forces in defense 
of United States' interests throughout the world mandate 

the development of fuel efficient, high speed, minimum man- 

b b 

ned ships for the Navy. The Captured Air Buddie (CAB) Sur- 
face Effect Ship (SES) holds great promise in providing an 
answer to these problems. 

A ship with minimum manning traveling at the high speeds 
envisioned for the SES would require highly skilled operator 
personnel. It would be necessary to provide these personnel 
with real time operational training similar to that which 
pilots receive in aircraft simulators. This would test the 
man in responding to failure situations that may be encoun- 
tered during actual operation and train personnel in control 
techniques that are unique to SES type vehicles. A Real Time 
Simulator could also be used as a means to test hardware used 
in SES craft and provide a design tool for the development 
of future modifications. 

A simplified, non-linear five degree of freedom Real Time 
Simulator (RTS5D) for the 3K-SES has been developed by T. S. 
Nelson [Ref. 1] such that a continuously observable real time 
solution is interfaced \vith man-generated control. Certain 
"man- in- the- loop” experiments have been conducted using the 
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RTS5D that have demonstrated limitations in that model. It 
was necessary to: 1) refine ship dynamics; 2) introduce sim- 
plified propulsion dynamics; 3) constrain thrust angle in 
order to remain within the bounds of programmed dynamic res- 
ponse (avoid broaching) ; 4) refine the software to provide 
faster computer iteration time; and 5) change the hardware 
to provide the operator with speed control as well as turning 
contro 1 . 

Accordingly, changes have been made in the RTS5D model 
to more accurately reflect the nonlinear forward drag 
characteristics of the SES. A speed controller has been 
designed to allow operation of the 3K-SES within a speed 
range of 40-60 knots and limits have been placed on the 
thruster deflection angle in order to prevent broaching, a 
condition which is not allowed under the existing dynamic 
model. Software refinements have been incorporated to im- 
prove computer iteration time. 
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II. EQUATIONS OF MOTION 



A. COORDINATE SYSTEMS AND ASSUMPTIONS 

The RTS5D model of the 3K-SES used simplified equations of 
motion. For completeness of comparison the development of 
these equations is repeated here. The rational for requiring 
simplification was to start with differential equations of 
motion that could be utilized effectively in a real time simu- 
lation. To satisfy this requirement a model was selected that 
would use point source of drag forces consisting of only a 
few lumped-coefficient terms. The equations used are based 
upon those developed in Ref. 2 iv r hich follow this concept. 
Additional assumptions used in the simplified 5 DOF equations 
were as follows: 

1) All accelerations are measured at the center of 
gravity. 

2) Vehicle is "free- to-heave" . 

3) All cross coupled moments of inertia are zero. 

4) All cross products of angular velocities in force 
equations are zero. 

5) All roll and pitch angles used in force equations 
are subject to small angle approximations. 

6) Calm water conditions. 

7) Mass and mass distribution of vehicle are constant. 

8) Effect of changes of aerodynamic force are neglected. 
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Utilizing the above assumptions, the simplified 5 DOF 
equations of motion are developed using the coordinate sys- 
tems shown in Fig. 1 and Fig. 2, the force vector diagrams 
shown in Fig. 3, Fig. 4A, Fig. 4B, and the following 
equations from Newton's laws of motion: 



SURGE 


m (u - 


vr) 


= X 


SWAY 


m(v - 


ur) 


= Y 


YAW 


! z r 




= N 


PITCH 






= M 


ROLL 


: x p 




= K 



where : 

m = mass of the rigid ship 

u = velocity in x-direction (SURGE) 

v = velocity in y-direction (SWAY) 

r = angular velocity about the z-axis 

p = angular velocity about the x-axis 

q = angular velocity about the y-axis 

= moment of inertia about the z-axis 

Iy = moment of inertia about the y-axis 

1^ = moment of inertia about the x-axis 

X = summation of forces in x-direction 

Y = summation of forces in y-direction 

N = summation of moments about the 
z-axis 

M = summation of moments about the y-axis 
K = summation of moments about the x-axis 



' slugs 

1 ft/sec 

' ft/sec 

' deg/sec 

' deg/sec 

' deg/sec 

' lb -ft 2 
m 

T lb -ft 2 
m 

' lb -ft 2 
m 




' ft- lb £ 

' ft-lb f 
' ft-lb £ 
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Additionally the following navigation relationships are 
utilized (see Fig. 1). 



X 

0 


= U COSlj) 


- v sin^ 


Y 

0 


= u sin<ji 


+ V COS\p 


2 

V 

s 


2 2 
= u + V 




B 


= tan 


(-V/U) 



where : 

ij> = heading angle 

6 = drift angle 

5 = thrust vector angle 

Also , 



X 

0 


= fX dt 
0 


u = 


/ u 


dt 


Y 

0 


= /Y dt 
0 


V = 


/V 


dt 


r 


= fr dt 


ip = 


fr 


dt 


P 


= / p dt 


* = 


fp 


dt 


q 


= fq dt 


9 = 


f q 


dt 
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Note: 1) 



ip = 0 when vehicle heads paralled to + X„. Tr 
axis. NAV 



2 ) 



Vehicle in figure above shown with 
thrust vector angle 6, positive ip, 
8, negative sway velocity v, i.e., 
turn 



negative 
negative 
in a right 



Figure 1 

Definition Of Coordinate System (Part I) 



\ 
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a 



r ^ 





Note : 



Direction for positive yelocity and 
angular rates are shown 



Figure 2 

Definition Of Coordinate System CPart II) 
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Note: 1) p and p' are equivalent point force centroids. 

2) Sj, 6 g , <5g, are thrust vector angles and 

are not required to be equal. 

3) Tj, Tg, Tg , are thrust magnitudes and are 

not required to be equal. 

4) Force directions shown are for a right turn 

Figure 3 

Surface Effect Ship CTop View) 
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+ z 




Note : 



1) Right Turn Forces Acting, i.e., Sj, S g , <5 g , 
6^g are negative 

2 ) T SIDE “ T 7 5ln5 y +T 8 SinS 8 tT 9 Sin6 9 +T 1 0 sinS 10 



Figure 4A 

Surface Effect Ship (Stern View) 
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+ z 




Note: TpQRW T^cos<Sy + TgCosog + TgCos6g + T-^QCOS(5^Q 

(see Fig. 3) 



Figure 4B 

Surface Effect Ship (Side Vievvr) 
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B. FORCES AND MOMENTS 



1- Surge Forces [see Fig. 4B) 

The forward acceleration was determined to be genera- 
ted by the summation of thrust vectored in the forward direc- 
tion T^ orw . The resultant acceleration of the SES was 
assumed to be opposed by a retarding force exhibiting a 
"velocity- squared” characteristic. Additionally, Newton's 
laws of motion specified a "v.r" product that contributed 
to the drag component of the equation when the vehicle was in 
a turn (contrifugal force reaction term) . The simplified 
equation of motion for the surge acceleration follows from 
a summation of forces in the u direction as depicted in 
Fig. 4B. 

u = (T forw /ra) ' ( C DX u2/m) + vr 

2. Sway Forces (see Fig. 4A) 

The side acceleration of the vehicle was analyzed 
to be the result of the summation of the thrust vectored 
parallel to the vehicle's y-axis and a retardation "velocity 
squared" phenomena such as that produced by the cross flow 
drag term of the sidewall as described in Ref. 3. This 
retardation force is augmented by a "u*r" product as speci- 
fied by Newton's Law. 

v ' T side /m • C DY v|v|/m ' ur 
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7 

Note the requirement to force the C^v /m term to 
maintain the sign of the sway velocity such that both left 
and right turns may be computed. 

3. Yaw Moments (see Fig. 3) 

The yaw acceleration was a turning moment summation 
which was defined to be influenced by three primary forces 
operating over moment arms. It is of interest to note that 
the third term is an "added mass" term which is considered 
to operate on the same pressure point as the Cpyv|v[ term. 
The term T ^ is a summation of moments generated by the 
four thrusters as shown in Fig. 3. 

T yaw = S 4 T 7 cos5 7 + S 3 T 8 cos 5 8 ‘ S 1 T 9 C0S 6 9 

S 2 T 10 C0S 5 10 “ £ o T 7 Sin 6 7 ' £ o T 8 Sin 6 S 
* £ o T 9 Sin 5 9 • *o T 10 Sin 5 10 

r = C DY Z w V > V I /I Z + A 22 UVZ w / ' I Z + T YAW /: Z 

4. Pitch Moments (see Fig. 4B) 

The pitch acceleration was assumed to be the summa- 
tion of moments generated by forward thrust, T£ Qrw , the 
buoyancy of the sidewalls of the SES, F^ and F- and a 
vertical force generated at the bow of the vehicle, F-. 

This vertical bow force was defined to be a lumped parameter 
term which modeled the reaction force due to plenum pressure 
acting against the bow seal. 
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z 



d 



A 

w 



1 




Z 



xl 



A 

w 



1 



Z 



*2 



w 

e 



A 



w 



2 




P 



cr 

o 



z 
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z 



3 



q 



- average draft of bow seal 

= £ x 5 'x 1 + CC* X tane)/2)t x 

= i i * ((a tane)/2t 

A A ^ A A 

= A Z pg 

W 1 X 2 

— A £ p g 

w 2 x 2 

C DpPb W e^d * *31 tan 6 ' ) 



= average draft of sidewall 
= plenum pressure 

= average wetted sidewall area of the stern 



= width of one sidewall 



= width of bow seal 

= average wetted sidewall area of the bow 

= added mass coefficient 
= water density 
= 32.3 ft/s 2 
= lever arm of bow seal 
= lever arm of buoyancy force 



(T. t + F.jU, + F~z 
v forw p 3 pi 2 o 




z 

z 



F l 2 3 ' A 33 Uq ' )/I Y 
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Note the added mass term A^uq which was required 
to provide damping to the pitch moment and is specified 
in Ref. 2. 



It was found that the response of the SES to a 

pitch perturbation without the added mass component was 

undamped and approximately sinusoidal. Unlike the yaw 

equation where the added mass term was a small contributor 

to damping, the pitch added mass term was found to be of 

significant importance in modeling the known pitch motion. 

5. Roll Moments (see Fig. 4A) 

The roll acceleration equation was analyzed to be a 

summation of moments generated by buoyancy forces of the 

port and starboard sidewalls, F and F , the thrust 
^ ’ sp ss’ 

vectored parallel to the y-axis of the vehicle, T s ^ e , a side 
force C^y v|v|, and a lumped coefficient vertical force, 

F . The vertical force F was defined as the force 



sw sw 

generated in a turn due to the sidewall curvature (dead 
rise angle) acting against the cross flow of \vater . 

F = pg A idp 

sp 6 wp 

F = pg A 2ds 

ss & ws 

F = v|v| 

sw DZP 1 

Z, = average draft of SES sidewall 



Z 

l 



dp 

ds 




l tan <f> 
sw r 

2, tan d> 
sw 
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= average wetted area port 
= average wetted area starboard 




P 




F )£ 
ss sw 



T . , £ 
side p 



+ C 



DY 



v v 



£ 

z 



- F £ 
sw sw 



A 31 up)/I x 



Again note the added mass term A,,up/I Y which was 
found to be essential in modeling the damping phenomena of 
the roll motion. 

In summary, the simplified 5 DOF equations of motion 
used in RTS5D for the 3K TON SES are: 



SURGE u = T forw /m - C DX u 2 /m + vr (Fig. 4B) 

SWAY v = T si d e / m ’ c D y v M/ m = ur (Fig- 4A) 

M* r = T yaw /I z - C DY v|v| £ w /I z + A 22 uv£ w /I z 

(Fig. 5) 

£U£H « ' (T forw t p * F 3 i 31 * F 2 *3 - C DX u 2 tz 

- F 1 £ 3 - A 34 uq)/I y (Fig. 4B) 

Mi P = (CF sp - F ss )* sw - T side z p * C BY v|vU 2 

- F sw l sw - A 33 Up)/I X (Fi «- 4A) 
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C. PARAMETER IDENTIFICATION 



These equations are an extension of the 3 DOF flat turn 

SES model developed by Gerba and Thaler in Ref. 4. The 

identification of craft parameters C..,,, C and Z is 

r DX DY w 

described in Ref. 4 and repeated here for completeness. 

The surge drag coefficient is determined by selecting 

a steady state turn condition, where u, v, r, and m 

are known. 



u = 0 = 



forw 

m 



C DX U 

m 



+ vr 



from which 



C 



DX 



forw 



u 



+ mvr 



The sway drag coefficient C^y is determined by the same 
steady state turn condition which requires 



v 



0 = 



side 

m 



"DY 

m 



v v 



ur 



where 



C 



DY 



T 



s ide 
v | v 



mur 



The sway drag moment arm follows utilizing the yaw accel- 
eration equation 



r = 0 



yaw + _DY 



C - VIV '\ + A 2 2 uv V 



I. 



I, 



I. 
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which yields 



-T 



w = 



yaw 



J DY 



v|v| + A 22 uv 



The addition of the pitch and roll equations introduces 

additional lumped coefficients C^p and C DZp . These are 

easily solved using known constants 1,1, C™, 1,1 

° w* p’ DY’ z* sw’ 

w p, , l, and steady state values of F , F , F, , F_, 

e ’ r b d } spss’12 

^forw’ T side’ u > v> anc * r ‘ s ig n ifi- cant t0 note that 

in a steady state condition the angular roll rate, p, and 
angular pitch rate, q, are both required to equal zero. 

This is in contrast to the yaw rate, r, where a finite 
steady state value is desired- in a turn. Thus the added 
mass terms of A^-uq and A^ p up in the pitch and roll equa- 
tions are not utilized in the determination of craft para- 
meters; their function is strictly confined to damping of 
their respective accelerations. Therefore, utilizing the 
pitch equation with known steady state turn values. 

q = 0 = ^ T forw £ p + C DP P b W e^d ' l 51 tan 0) *31 

+ F 2 £ 3 ' C DX U 12 " F 1 Z, 3' )//I Y 

yields 

_ C DX u2jl z + F l £ 3 ' F 2 £ 3 ' T forw *p 
DP p bV £ d - *31 tan 0J '*31 
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The roll equation under steady state turn conditions 
yields 



P ’ 0 = ((F sp ' F ss )i w ■ hide “p + C DY v|v l iz 



C DZP v l v l *5W )r z 



from which 



"DZP 



■ C(F sp - F „J*w T side l p * C DY hs'hs'h 57 



where 



(v |v | £ ) 

v ss 1 ss 1 sw ; 



F sp 


— 


port sidewall buoyancy in pounds 


(Fig. 


4A) 


F 

ss 


= 


starboard sidewall buoyancy in pounds 


(Fig. 


4A) 


T 

1 side 


= 


thrust vectored parallel to vehicles y 
in pounds 


-axis 

(Fig. 


3) 


V 

ss 




sway steady state velocity in feet/sec 






C DY 


= 


sway drag coefficient for steady state 


condition 


l 

w 


= 


sway drag moment in feet 


(Fig. 


3) 


l 

p 


= 


effector thrust moment arm in feet 


(Fig. 


4A) 


l 

z 


= 


surge drag moment arm in feet 


(Fig. 


4A) 


l 

sw 


= 


sidewall buoyancy force moment arm 
in feet 


(Fig. 


4A) 
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III. RTS5D VALIDATION 



A. DBSIMSD BENCHMARK 

Validation of the RTS5D was attempted "using the data 
base program (DBSIMSD) [Ref. 5] as a benchmark. However, 
new guidance in the proper use of the data base program con- 
tained in Ref. 6 revealed the following errors in the use of 
the DBSIMSD. 1) Initial conditions had been improperly set 
and 2) Data used for validation for the 56 knot run at 
effector deflection angles greater than 10° was not correct 
because of the thrust inlet broach conditions which existed 
using the DBSIM5D. These conditions were not accounted for 
in the RTS5D model. The noted discrepancies in the opera- 
tion of the DBSIM5D necessitated a revalidation of the 
RTS5D. 

B. VALIDATION RESULTS 

Using the information contained in Ref. 6, the DBSIM5D 
program was again used as a benchmark to validate the RTS5D 
simulation. A DBSIMSD sequence of runs with initial forward 
velocities equal to 40 knots, 50 knots and 60 knots with 
step effector angles of 5° and 10° (at 40 knots), 5°, 10° 
and 15° (at 50 knots), and 5° (at 60 knots) were compared 
to identical runs of the RTS5D . Thrust effector angles 
were kept below a maximum that would have caused a broach 
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condition. Broach conditions are explained in Chapter IV. 

The first peak overshoot value and "Quasi Steady State" 
values for the variables u (surge) , v (sway) , r (yaw rate) , 

4> (roll) , and 0 (pitch) were used as a basis for comparison. 
The results are shown in Tables I - 1 1 1 . Significant error 
in forward velocity occurred in the 50 knot test. Also, it 
was noted that in 40 knot test with 5° thruster deflection 
angle that the DBSIM5D actually reached a "Quasi Steady 
State" forward velocity greater than the initial value upon 
entering the turn. These conditions led to the development 
of the model of the RTS5D with negative drag characteristics 
in sea state. Differences in the other measured parameters 
were addressed by parameter adjustment on the new model in 
order to more closely match values obtained with the DBSIM5D. 
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TABLE I 



RTS5D and DBSIM5D Performance Test at 40 Knots 



effector 

angle 


1st pk 


DBSIM5D 


in 

Qss 


RTS5D 

degrees <t> 

1st pk 


Qss 


5° 


1.01 




. 742 


.17 




.12 


10° 


1.46 




1.39 


.31 




. 24 




1st pk 


$ 


in 

Qss 


degrees 

1st pk 


* 


Qss 


5° 


1.44 




1.40 


1.27 




1.21 


10° 


1.48 




1.48 


1.28 




1.20 




1st pk 


u 


in 

Qss 


ft/sec 

1st pk 


u 


Qss 


5° 


67.4 




70.9 


n/a 




67.03 


10° 


n/ a 




66.3 


n/ a 




65.98 




1st pk 


V 


in 

Qss 


ft/ sec 

1st pk 


V 


Qss 


5° 


1.42 




1.37 


2.84 




2 .42 


10° 


3. 20 




3.26 


3.83 




3.41 




1st pk 


r 


in 

Qss 


degrees/sec 
1st pk 


r 


Qss 


5° 


.396 




.403 


. 58 




. 33 


10° 


.953 




.959 


.98 




.66 
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TABLE II 





RTS5D and DBSIM5D 


Performance 


Test at 50 


Knots 


effector 

angle 


DBSIM5D 

1st pk 


in degrees 
Qss 


RTS5D 

* 

1st pk 


Qss 


5° 


.424 


.446 


. 27 


.19 


10° 


.757 


.886 


.49 


.38 


15° 


1.02 


1.4 


. 72 


. 57 




e 

1st pk 


in degrees 
Qss 


6 

1st pk 


Qss 


5° 


.862 


. 892 


1.16 


1.15 


10° 


. 943 


1.13 


1.15 


1.16 


15° 


1.12 


1.35 


1.22 


1.17 




u 

1st pk 


in ft/sec 
Qss 


u 

1st pk 


Qss 


5° 


n/ a 


83.6 


n/a 


83.8 


10° 


n/a 


78.7 


n/a 


82.51 


15° 


n/a 


69.2 


n/a 


80.58 




V 

1st pk 


in ft/sec 
Qss 


V 

1st pk 


Qss 


5° 


1.31 


1.30 


3.55 


3.02 


O 

o 


2.72 


2.70 


4.79 


4.26 


15° 


4.05 


4.33 


5.69 


5.207 




r 

1st pk 


in degrees/sec r 

Qss 1st pk 


Qss 


5° 


.368 


.365 


.73 


.41 


h- * 
O 
O 


. 864 


. 870 


1.22 


.82 


15° 


1.35 


1.44 


1.66 


1.26 
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TABLE III 





RTS5D and 


DBSIM5D 


Performance Test at 50 Knots 




effector 

angle 


DBSIM5D 




RTS5D 










in degrees 


♦ 




3° 


1st pk 
.141 




Qss 
. 200 


1st pk 
.39 


Qss 
. 274 






9 


in degrees 


0 




5° 


1st pk 
. 521 




Qss 
. 583 


1st pk 
1.14 


Qss 

1.07 






U 


in ft/sec 


u 




5° 


1st pk 
n/a 




Qss 

97.8 


1st pk 
n/a 


Qss 

100.29 






V 


in ft/sec 


V 




5° 


1st pk 
2.03 




Qss 

1.84 


1st pk 
4.26 


Qss 

3.62 






r 


in degrees/sec 


r 




5° 


1st pk 
. 532 




Qss 

.494 


1st pk 
.87 


Qss 

.487 
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IV. MODEL DEVELOPMENT 



A. DRAG CHARACTERISTICS 

1 . Surge Equation Model 

The RTS5D model assumed a velocity squared drag re- 
tarding characteristic for the 3K-SES which led to the equa- 
tion for surge acceleration that was presented in Chapter II. 
The drag characteristics from Ref. 7 for the 3K-SES that is 
simulated by the DBSIM5D are shown in Fig. 5. The data base 
program was run using parameters that caused it to exhibit 
a retarding force corresponding to the curve at sea state 
four, therefore this curve was used to remodel the RTS5D. 

The new simplified equation of motion for surge acceleration 
became 

u = ( T f 0rw /m') - (FSD/m) + v . r 

where FSD is full scale drag which is defined by the curve 

in Fig. 5. T is the sumation of the four thrusts assumed 

rORW 

in this model to be applied at the same point along the 
center line of the craft. 

2 . Implementation 

FSD was implemented into the RTS5D model by dividing 
the curve into four separate regions and defining FSD within 
each region. In region I the retarding force was approxi- 
mated by a normalized velocity cubed term. In regions II 



36 




I 

o 

bjD ^ 

cti X 
U CO 

a & 



0 ) 

3 

W) 

•H 



37 



Full Scale Drag Curves 



and III straight line approximations were used. In region 
IV a polynomial curve fit was used. The equations defined for 
each region are listed below. 

Region I - ( 0-18 knots) 

FSD = C M1 ( u/ u max )3 
Region II - (18-22 knots) 

FSD = CONSTANT^ - u 

Region III - (22-36 knots) 

FSD = C DX2 u - CONSTANT^ 

Region IV - (above 3 6 knots) 

FSD = C DX ^ u 2 + C Dx j^/ u 1 * 5 

3 . Identification of Surge Drag Coefficients 

In region I, is equal to the first peak value 

of drag. In regions II and III, C^,. and Cp X? are equal to 
the slopes of the drag curve and the constants were solved 
for by applying the boundary conditions at the intersections 
of the adjacent regions. The values of and were 

computed by taking F.SD at two different known values of u 
and solving for C DX3 and C~ x4 simultaneously. 

B. LINEARIZED MODEL ANALYSIS 

Since the surge model of the RTS5D is extremely non- 
linear over the complete speed range, analysis is very com- 
plex. Therefore the model was linearized and it’s performance 
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analyzed at various points on the FSD curve. Linearization 
of the model was accomplished by computing the linear 
coefficients through partial differention with respect to 
surge in each of the four region equations and evaluating 
these terms at various operating point surge velocities (u ) . 
The linearized model is shown in Fig. 6 and the linear coef- 
ficients for the surge drag force FSD are shown below. 

Region 

3 C DX1 ( u 0 / u max )2 
" C DX5 
°DX2 

2 C DX3 U o ' 1,5 C DX4/ u o 2 ’ 5 

A map of the linearized surge directional roots are 
shown in Fig. 7. The characteristic roots for each operating 
region are shown coded Csee legend) . Note that only Region I 
and Region IV roots are u q dependent and therefore move as u q 
changes. The direction of change along the real axis is 
noted by the arrows shown. Region II and III roots are fixed 
in locations and independent of u . 



/ 



I 


FSD (u 0 ) ■ 




/ 


II 


FS X> ■ 




/ 


III 


PSD (u n ) = 




0 


IV 


FSD, \ = 
(u o> 
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Figure 6 

Linearized Surge Model 

Region IV of Fig. 7 is of particular interest since the 
ability to operate in this speed range is a primary consi- 
deration in the development of the SES. It can be seen 
from the map of the roots that within this region, at speeds 
from 36 to almost 60 knots, poles exist in the RHP. This 
is a result of the negative drag slope where the retarding 
force actually decreases with an increase in speed, the 
magnitude of variation being dependent on sea state. Region 
II also exhibits this negative drag characteristic but since 
it is of such a narrow speed range, operation in this re- 
gion was not considered in this design. Further modeling 
changes for the other degrees of freedom were based on 
Region IV response. In Regions I, III and in Region IV above 
60 knots, the system will be stable at the speed u q where FSD 
is equal to commanded thrust. 
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Figure 7 

Linearized Surge Root Map 



If it were necessary to operate at a speed within the 
negative drag slope characteristic, it would be necessary to 
actively modulate the thrust command in order to maintain 
a nearly constant speed, but this would require a great deal 
of operator attention since the negative drag slope in Region 
IV is one' of slowly divergent instability. For this reason 
a closed loop speed controller has been included in the data 
base model and therefore was also added to the RTS-51) model. 

C. SIMPLIFIED PROPULSION" DYNAMICS 

Before adding a speed autopilot to the model one addi- 
tional refinement was included in the design. Thrust command 
changes in the RTSSD as reported in Ref. I were instantan- 
eously entered into the equations of motion as thrust 
changes. This was not the case for the data base simulation 
program. Specific thrust levels are generated by specific 
water flow rates through water jet nozizles. These flow 
rates are nearly proportional to the pump speed that pro- 
duces them. Since the pumps are driven directly by gas 
turbine engines, a change in thrust occurs with a change in 
turbine speed. A thrust command translates to a power tur- 
bine speed command which as an output has power turbine 
speed achieved. This output is translated into thrust 
achieved. It was specified in Ref. 7 that one of the goals 
in the base line design was to have power turbine speed 
control response characteristics which are relatively inde- 
pendent of the power setting. Accordingly, in the design 
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of t.he auto pilot of Ref. 7 , t he simplified propulsion 
fa lock shown in Fig. 8 was used. The magnitude of the pro- 
pulsion system time constant, t was not specified. For 
the purpose of this analysis t was chosen to be 1 second. 



$ Fuel 



Command 


1 


<f T 




l +Ts 





Figure 8 

Simplified Propulsion System 



D. ANALYSIS OF LINEARIZED MODEL WITH SPEED CONTROLLER 

The complete linearized model with speed controller that 
is included in the model for the purposes of this study is 
shown in Fig. 9. The saturation effects of fuel control are 
ign.ored and the linear representation is shown as which 
is assumed to be unity. The controller includes propor- 
tional plus integral control. The proportional control was 
included to keep the ship speed error low while the integral 
control was included in order to zzero out long term control 
error . 
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Propulsion System and Speed Controller 



The characteristic equation for the system in Fig. 9 is 



S 3 + (1 + p SD( u )/ m ) S 2 + l/m (FSd' (u j + K) S + K/ =0 



Using the Routh Hurwitz Criterion for stability 




1 



1/m (FSD' (u } + K) 



S 



2 



( 1 + FSd' (u j/ m ) 



K/ m T 



S 1 1/m (1 + FSD^ u j/m) (FSd' (u j + K) - K/m r 



0 



S 



0 



K/ m r 



0 



For Stability K must be greater than 0 and greater than 



The values of controller parameters K and t can be initialized 
by the operator for either loosely controlled speed or tightly 
controlled speed for such situations as station keeping, 

UNREPS, or test operations. For comparison to the data base, 
these parameters were selected to achieve closest agreement 
between the values of surge achieved by the RTS5D and DBSIM5D 
models after the completion of a 360° turn with 50 knot initial 
surge veloctiy and a thrust effector angle of 15°. This test 



FSD '(u o ) r + »™a' ( , 
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condition for the nonlinear models was chosen for two reasons 
First, given fuel conservation as a necessity, 50 knots is a 
likely high speed initial condition because it is close to 
minimum drag. Second, this initial condition provided an 
opportunity to test the operation of the speed controller in 
the negative drag region of FSD since the decrease in speed 
going into the turn with u q =50 knots caused increased drag. 
The values of controller parameters that caused closest agree 
ment between the two nonlinear models, when tested in the 
linearized model showed instability in the linear sense. 

This is desirable in a turning mode for the nonlinear model 
because it allows speed to decrease and the turn to be com- 
pleted more quickly. Stiff speed control could be achieved 
by increasing the controller gains, but would normally be 
used only in straight ahead runs. 

Coupling of the equations of motion developed in Chapter 
II necessitated a repetition of the parameter identification 
procedure resulting in new values for C^y, £w, C^p and C^ 7 p. 
These values were adjusted for the best curve fit compared 
to the data base using the linearized model for the RTS5D 
in a 360° turn using 15° thrust effector angle with initial 
surge velocity 50 knots. An identical series of tests to 
those in Chapter III were conducted comparing the new model 
of the RTS5D to the DBSIM5D. Complete validation results 
are tabulated in Chapter VI. 
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E. EFFECTIVE THRUST EFFECTOR ANGLE 



As a result of the negative drag characteristics imposed 
on the model , RTS5D yaw rate (r) was high in the low speed 
test and low in the high speed test. This opposite error 
at both ends of the test speed spectrum was reduced by 
introducing a correction factor to the thrust effector angle 
proportional to yaw rate 

Z .. = Z - SLIP • r 
ef f 

and adjusting £w to increase r. This had the desired effect 
of reducing error at both ends of the test spectrum since 
the reduction of the thrust effector angle was greater at 
the high end. 

F. BROACH CONDITION FLAG 

As stated in Chapter I, broach response characteristics 
are not included in this model, therefore a broach condition 
warning flag was included in the design in order to alert 
the operator who could then take corrective action to avoid 
broaching . 

Broaching is caused by a water intake being exposed to 
air and is a function of surge speed, roll angle, pitch 
angle, drift angle and sea state. For the data base 3K-SES 
with broach boundaries defined at 101 air by volume, the 
broaching boundaries in calm water in terms of roll angle, 
pitch angle and drift angle are shown in Fig. 11 for 40 knots 
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and Fig. 12 for 60 knots. Reference 6 listed maximum thrust 
effector angles as a function of surge speed and air flow 
in order to avoid broaching. This condition is shown plotted 
in Fig. 10 for fixed plenum air flow rate. 



EFFECTOR 

ANGLE 




Figure 10 

Maximum Thrust Effector Angle As A Function 
Of Surge Velocity In Order To Avoid Broaching 



As stated in Chapter III the model developed in Ref. 1 did 
not flag the broaching condition. To implement the broach 
warning flag and display the other new parameters for the 
modified RTS5D the pilot graphic display described in Ref. 1 
was changed to incorporate the necessary data. Figure 13A 
depicts typical information displayed to the pilot under 
normal operating conditions. The new information to the 
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Inlet Broaching Boundaries at 40 Knots 



operator shown in the lower right of Fig. 13A is sea state 
(SEAST) , computer assist (CASST, on=l off=0) , controller 
gains (K,KK) and speed commanded (SPCOM) . if the operator 
exceeds the thrust effector limits (ZMAX) to prevent 
broaching (see Fig. 10), a flashing warning is displayed 
in place of the bottom line of information on the screen. 
This is shown in Fig. 13B. 
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Inlet Broaching Boundaries at 60 Knots 



VARIABLE DYNAMICS 



U (HEAVY) p*. 
UDOT (LIGHT) 



V (HEAVY) 
VDOT( LIGHT) 






R (HEAVY) 
RDOT (LIGHT) 






P (HEAVY) 
PDOT( LIGHT) 






Q (HEAVY) 
QDOT (LIGHT) 






SEA TRACK 





CONTROL INPUTS TIME NAVIGATIONAL DATA 



THRUST 


RUDDER 


DELT 


ELAPSED 


PITCH HEADING/RATE DRIFT 


ROLL SPEED 


300000 


-.20 


.084 


22.2 


.52 325/2.1 -.28 


.8 50.1 


T? 


T8 


T9 


T10 


SEAST CASST K KK 


SPCOM 


75000 


75000 


75000 


750000 4 1 2900 3 


50 



Figure 13A 

Pilot Graphic Display (Normal) 



VARIABLE DYNAMICS 



SEA TRACK 



U (HEAVY) 
UDOT( LIGHT) 



V (HEAVY) 
VDOT (LIGHT) 




R (HEAVY) 
RDOT( LIGHT) 






P (HEAVY) 
PDOT (LIGHT) 



Q (HEAVY) 
QDOT( LIGHT) 





CONTROL INPUTS TIME NAVIGATIONAL DATA 

THRUST RUDDER DELT ELAPSED PITCH HEADING/RATE DRIFT ROLL SPEED 

3000000 -.52 .084 22.2 .52 325/2.1 -.28 .8 50 . 1 



T7 T8 

~ B R 

/ 



T9 T10 SEAST CASST K 

0 A C H — 



KX 



SPCOM 



Figure 13B 

Pilot Graphic Display (Broach Condition Exists) 
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V. IMPLEMENTATION 



A. INTRODUCTION 

Given the existence of a Real Time Simulation of a SES , 
a decision was made to refine that model using more accurate 
equations of motion, improved computer iteration time and 
operator hardware modifications. A complete description of 
the RTS5D is contained in Ref. 1, therefore only changes 
are discussed here, and included are the introduction of 
negative drag in a sea state, propulsion dynamics, speed 
control and a broach condition warning in software. Hard- 
ware modification included changing the thruster control 
box to make it more accessable to an operator. 

B. REQUIREMENTS 

The following criteria were required of the remodeled 
RTS5D. 

1) The refined simplified equations of motion would be 
solved and the results output on a real t-ime basis. 

2) The solution would be subject to real time control 
efforts generated by an operator observing the output. 

3) Computer iteration, time was required to be less t-han 
Id <J ms . 

A more complex discussion of condition 3) is presented in 
Appendix A. 
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C . HARDWARE DESCRIPTION 



The hardware arrangement, in block diagram form is shown 
in Fig. 14. The only change to the RTS5 D was the redesign 
and relocation of the thruster console to the pilots seat 
where it can be more conveniently utilized by the operator. 
The thruster console uses "linear movement" potentiometers 
in place of previously used "dial" potentiometers in order 
to more realistically simulate thrust controls for the SES. 
The thruster console is shown iir. Fig. 15. 

D. SOFTWARE DESCRIPTION 

The remodeled software package was a FORTRAN IV digital 
program which consisted of a main graphics program and three 
major subroutines. It was necessary to make a third sub- 
routine out of the equations of motion computation that had 
previously resided in the main program prior to making any 
changes to these equations. It was discovered when the 
first program modifications were attempted that the compiler 
for the XDS-9500 computer was operating at its limit with 
the program as it then existed, therefore, before changes to 
the equations could be implemented, the program structure 
had to be modified. The failure mode was eliminated from 
the new model because of initial assumptions that the thrusts 
would be all applied at the same point.. Thrust failures 
should not be considered until the model includes broaching 
dynamics which is a major source of thrust loss. The flow 
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corn- 



chart for the remodeled RTS5D is shown in Fig. Id. The 
plete multiplexing algorithm is shown in a flow chart in 
Fig. 17. 
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Figure 14 

RTS5D Mod Block Diagram 
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A - Thruster T? 

B - Thruster T8 
C - Thruster T9 
D - Thruster TlO 
E - Speed Command Control 
F - Run, Hold, Restart 

Figure 15 
Thruster Console 



58 




Figure 16 



RTS5D Mod Program Flow Chart 
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DESCRIPTION OF INFO(I) 



INFO(l) u DYNAMICS 

INFO (2 ) u DYNAMICS 

INFO(3 ) v DYNAMICS 

INFO (4) V DYNAMICS 

INFO(5) r DYNAMICS 

INFO (6) r DYNAMICS 

INFO (7) p DYNAMICS 

INFO (8) p DYNAMICS 

INFO(9) q DYNAMICS 

INFO (10) q DYNAMICS 

INFO(ll) POSITION MAX VALUES 

INF0(12) POSITION MIN VALUES 



INFO ( 13 ) 
INF0(l4) 
INF0(15) 
INF0(l6) 
INFO (17 ) 
INF0(18) 
INF0(19) 
INFO (20 ) 
INFO (21 ) 
INF0(22) 
INF0(23) 



VELOCITY MAX VALUES 
VELOCITY MIN VALUES 
ACCELERATION MAX VALUES 
ACCELERATION MIN VALUES 
NAVIGATION DATA 
EFFECTOR/THRUSTER DATA 
SEA TRACK NAV PLOT 
SEA TRACK NAV PLOT 
PRESENT POSITION VARIABLES 
PRESENT VELOCITY VAREABLES 
PRESENT ACCELERATION VAR. 



Figure 17 

RTS5D MOD DISPLAY MULTIPLEX ALGORITHM 
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VI. RTS5D MODS I AND II RESPONSE CHARACTERISTICS 



RTS5D MOD I includes negative drag effects in a sea state, 
simplified propulsion dynamics, speed control and a broach 
condition warning flag. RTS5D MOD II is identical to MOD I 
except that it includes the correction factor for thrust 
effector angle developed in Chapter IV, Section E. 

The validation of the modified RTS5D was accomplished 
in the same manner as the RTS5D. CSee Chapter III) .. The 
results are displayed in the following tables. A complete 
360° turn comparison of the DBSIM5D, the RTS5 D and the modi- 
fied RTSS D at 50 knots, 15° thruster angle in shown in Fig. 

18. Additionally, a comparison of transient response was 
undertaken and the time of first peak for the variable v, r, 

$ and 9 for each model was graphed in Figs. 19-21 for speeds 
of 4Q , 5Q and 60 knots. Transient response characteristics 
are an important consideration in the development of a real 
time modeL and are discussed in Appendix A. The transient 
response of the RTSSD models was faster than that of the 
data base model with the exception of roll angle (4) . This 
impacted on the validity of a real time solution of the equa- 
tions of motion and the display of the results. 
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TABLE IV 



RTS5D Mod I 


and Mod II 


Performance 


Test at 10 


Knots 


effector 

angle 




RTSSD Mod I 


RTS5D Mod 


II 






<f> 


in degrees 






5° 


1st pk 
.52 




Qss 

.45 


1st pk 
.40 


Qss 
. 38 


O 

o 
1 — 1 


. 96 




.92 


.79 


. 76 






0 


in degrees 


0 




5° 


1st pk 
1.21 




Qss 

1.17 


1st pk 
1.20 


Qss 

1.18 


10° 


1.20 




1.18 


1.20 


1.18 






u 


in ft/sec 


u 




5° 


1st pk 
62.56 




Qss 

62.64 


1st pk 
62 . 58 


Qss 

67.65 


O 

o 


n/a 




61.37 


n/a 


61.71 






V 


in ft/sec 


V 




5° 


1st pk 
2.62 




Qss 

2.47 


1st pk 
2.33 


Qss 

2.27 


10° 


3.53 




3.52 


3.23 


3.20 






r 


in degrees/sec r 




5° 


1st pk 
.69 




Qss 
. 59 


1st pk 
.57 


Qss 
. 54 


10° 


1.18 




1 . 22 


1.04 


1.09 
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TABLE V 



RTS5D Hod I and Mod II Performance Test at 50 Knots 

effector 

angle RTS5D Mod I RTS5D Mod II 







<f> 


in 


degrees 


* 




5° 


1st pk 
. 50 




Qss 

.39 


1st pk 
.43 




Qss 
. 39 


10° 


.90 




.82 


. 83 




. 78 


15° 


1.31 




1.36 


1.24 




1.19 






6 


in 


degrees 


0 




5° 


1st pk 
1.23 




Qss 

1.19 


1st pk 
1.23 




Qss 

1.20 


10° 


1.24 




1.19 


1.23 




1.20 


15° 


1.24 




1.18 


1.23 




1.19 






u 


in 


ft/sec 


u 




5° 


1st pk 
n/ a 




Qss 

82.73 


1st pk 
n/ a 




Qss 

82.98 


10° 


n/ a 




78.68 


n/a 




79.55 


15° 


n/ a 




68.30 


n/ a 




72.75 






V 


in 


ft/sec 


V 




5° 


1st pk 
2. 57 




Qss 

2.30 


1st pk 
2 . 39 




Qss 

2.29 


10° 


3.48 




3.33 


3.25 




3.25 


15° 


4.13 




4.29 


3.98 




4.00 






r 


in 


degrees/ sec 


r 




5° 


1st pk 
. 57 




Qss 

.39 


1st pk 
. 49 




Qss 

.41 


10° 


.96 




.85 


. 88 




. 86 


15° 


1.31 




1.62 


1.25 




1.41 



63 



TABLE VI 



RTS5D Mod I and Mod II Performance Test at 60 Knots 



effector 



angle 




RTS5D 


Mod 1 


RTS5D 


Mod 


ii 






in 


degrees 


* 




5° 


1st pk 
. 53 




Qss 
. 39 


1st pk 
.46 




Qss 

.41 


5° 


1st pk 
1.23 


0 


in 

Qss 
1 . 20 


degrees 

1st pk 
1.24 


9 


Qss 

1.20 


5° 


1st pk 
n/a 


u 


Qss 

100.21 


1st pk 
n/a 


U 


Qss 

99.59 


5° 


1st pk 
2.64 


V 


Qss 

2.29 


1st pk 
2.49 


V 


Qss 

2.35 


5° 


1st pk 
.52 


r 


Qss 

.32 


1st pk 
.46 


r 


Qss 
. 36 
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z = 15° (effector angle) 



10 second intervals 

O RTS 5 D 
A DBSIM 5 D 
- RTS 5 D MOD I 
* RTS 5 D MOD II 

Figure 18 

360 Turn Comparison of Four Models 
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TIME 
1st pk 
(SEC) 



11 “ 
10 " 

9 “ 

8 



* 



A 

o 



X 



7 - 

6 .> 



A 

o 



5 



X 



X 



4 - 



3 " 



2 



1 .. 



o 

A 



% • ■■ ■ ■ ■ » 
v r 0 

u q = 40 knots 

Effector angle = 5 



«- 

9 



0-= RTS5D 

RTS5D Mod II 
X= DBSIM5D 



Figure 19 

Response Time Comparison at 40 Knots 
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st 

(SE 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 



X 



X 



A 



A 



o 



o 



X 



A 



X 



o 



A 



■ « i t | 

v r 0 9 

u q = 50 knots o= RTS5D 

Effector angle = 5° Mod II 

X = DBSIM5D 

Figure 20 

Response Time Comparison at 50 Knots 



13 

12 

11 

10 

9 

3 

7 

6 

5 

4 

3 

2 

1 



X 






X 



A 



o 



A 



o 



A 



0 



X 



V 



A 



r 0 0- 



u = 60 knots 
o 

Effector angle = 5 

Figure 21 

Response Time Comparison at 60 Knots 



o = RTS5D 
A = RTS5D Mod II 
X = DBSIM5D 



68 



VII. CONCLUSIONS 



The RTS5D has been shown in Section VI to be a viable 
input/output model of a 3K-SES by closely matching its output 
characteristics for a given input to those of the data base 
mode 1 . 

The refinement of the right side of the simplified equa- 
tions of motion, especially the inclusion of the negative 
drag characteristic in a sea state has made the model more 
accurate over a wide speed range and gives the operator, upon 
program initialization, the ability to choose different drag 
characteristics caused by sea state. 

The model was tested in a negative drag region with 
automatic speed control and simplified propulsion dynamics 
and it was demonstrated that by properly adjusting controller 
gains that the output characteristics of the DBSIM30 were 
closely matched. 

Previously conducted thrust failure analysis with man 
generated corrective action has been shown to be invalid 
until broach dynamics are included into the model. Previous 
tests allowed maximum thrust effector angle correction regard- 
less of this consideration. The modified version of the RTS5 D 
does not include complete broach dynamics but does warn the 
operator, if he should exceed a maximum thrust effector angle 
for a given forward speed. 
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The 3K-SES may be operated in the negative drag slope 
speed region approximately C36-S8) knots without the aid of 
a speed controller but would require constant manual modula- 
tion of the thrusters. Since minimum manning is one of the 
criteria in. 3K-SES design, a speed controller is necessary 
to free operator attention to other duties when operating 
in this speed range. 

Anotfrer consideration of the implementation of the non- 
linear drag curves is that for a given sea state, two values 
of minimum drag can be determined, one in Region IV and one 
between Region II and Region III. Knowledge that drag can 
be lower in a particular operating region for fuel conserva- 
tion is of prime importance. 

It is recognized that a time delay exists in the solution 
of the equations of motion and a time delay exists in the 
graphical display of the solutions to the operator. Real 
time interaction between the man in the loop and the computer 
is achieved by setting the integration step size equal to 
th.e computer loop iteration time. Real, time considerations 
are further discussed in Appendix A. 
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VIII. RECOMMENDATIONS 



The following recommendations are made concerning future 
improvements in the RTS5D model. 

1. Further refinement of the right hand side of the 
equations of motion is necessary in order to more closely 
match the characteristics of the data base model. 

2 . Implement software broach dynamics in the model with 
thrust applied differentially at different lever arms Csee 
Fig. 3) so that thrust failures and effectiveness of correc- 
tive action may be evaluated. 

3. Rewrite the program such that the AGT-10's would 
provide asyncronous parallel graphics ("existing equipment) 
or implement the model on a more modern system which uses 
parallel processing and improved graphics. 

4. Introduce the sixth degree of freedom Cheave) as a 
parallel operation. Since heave time constant is much faster 
than the rest of th.e system dynamics, real time heave analysis 
could be conducted simultaneously if cross coupling is neg- 
lected. 

5 . The system is designed for man in the loop interface, 
so further operator te. sting is recommended as a valuable 
source of information for future refinement of the model. 
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APPENDIX A 



REAL TIME ANALYSIS 

T. S. Nelson, in Ref. 1, accomplished real time solution 
of SES equations of motion with the use of a 1000 Hz clock 
to measure computer iteration time and by using that time as 
an integration step size for rectangular integration to solve 
the simplified dynamic equations. He concluded, by comparing 
th-e solutions of these equations to that obtained by using a 
Runga Kutta variable integration step size, that the computer 
iteration loop time must stay below 100 ms in order to main- 
tain at least 3% accuracy of solution. He further developed 
a multiplexing algorithm in order to output these solutions 
to an operator in the form of graphics on a real time basis. 
This model can be invalidated as a real time model in two 
ways. First, in order to have real time, the integration 
step size used must be equal to the computer iteration time 
and in this model, because of the multiplexing algorithm, 
each loop time is different. Tfie integration step size, 
dt determined by measuring loop N is used to solve the equa- 
tions in loop N+l, which has a different loop time, there- 
fore the solutions are not exactly real time. Secondly, 
the multiplexing scheme makes the model into a sampled data 
system with sample time (Tsl approximately 2 seconds. Even 
if the first condition did not apply and real time solutions 
to the equations could be assumed, then the output to the 
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operator would be sampled real time data. Given the fastest 
time constant for the RTS5D simplified equations (approxi- 
mately 1.3 seconds for pitch rotational motion), it can be 
seen that the Nyquist Criteria is clearly violated. Yet 
even with these obvious limitations, it has been demonstrated 
that the RTS5D can closely reproduce the output characteris- 
tics of the DBSIM5D and therefore is a useful model. Surge 
speed as well as sway, yaw and roll (Figs. 21-23) have much 
longer effective time constants. It is recognized that 
there is a delay between the real time of any event and the 
observation of that event, therefore "Real Time" as used in 
this report is a matter of degree, or more specifically 
how much delay time between solution and display can be 
tolerated. Given the scale of graphics presented to the 
operator in this model and the speed with which the SES 
responds, the RTSS D is an acceptable model, however it 
would be desirable to reduce the sample time T^ and make the 
integration step size equal to the computer loop iteration 
time for each iteration. 

An algorithm to set the integration step size equal to 
the computer iteration, time of the current loop was imple- 
mented as a software design, change to the RTSS D model. This 
was accomplidred by defining the integration step size 
as an array DELTCN) with N equaL to th.e total number of itera- 
tion, s ne-cessary -t o display all of the graphical information 
to t'he operator. Each loop iteration time was ire asu red and 
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stored in. this array during th-e first pass through, the 
mu l"t iplexing algorithm. DELTCMJ was initialized at z.ero 
so that time would not increment and integration would not 
be performed during the first pass through the multiplexing 
algorithm. On subsequent passes through the multiplexing 
algorithm the integration step size was equal to the loop 
time in which it was used. 

This timing algorithm when tested in the RTS5 D model 
was found to accurately represent real time but the devia- 
tion from previously used integration step size values 
caused instability in both the pitch and roll dynamics. 

This needs further investigation. 

Sturgeon, in Ref. 8, demonstrated the value- of using 
parallel array processors to reduce computational time in 
real time simulations. Such a model could be designed for 
the RTSSD using the existing equipment. The AGT-10's which 
are- used for displaying the output can interface with the 
XDS-93O0 main computer in two ways. The multiplexing scheme 
of the RTS5D uses a "hand shaking" method with a GRAPHO or 
TEXTO subroutine call in each loop. It is precisely the 
inclusion of graphics in each loop that causes problems for 
real time presentation. The solution of the equations requires 
only approximately 35 ms. Each GRAPHO or TEXTO call, with 
its associated "hand shaking" time, makes up the rest of the 
loop time, approximately 40 ms more. It can be seen that by 
introducing two such calls in one loop that the 100 ms time 
restriction would be exceeded. 
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An improvement to this method would be to use the second 
interface between the AGT-lQ's and the XDS-9300 which is 
that the AGT-10's have a direct connection to some memory 
locations in XDS-93 00. By solving the dynamics for th.e SES 
and storing the data in these memory addresses in the form 
of arrays, the AGT-lCT’s could be made to operate asyncronously 
in parallel as array processors by using a process known as 
"cycle stealing" to fetch the data from the XDS-9300 and to 
display this on the screen for the operator. This process 
would have the advantage of reducing and making more uniform 
the computer iteration time used as the step size for the 
rectangular integration in the solution of the SES equations 
of motion. The sample time would also be reduced because all 
of the time wasted in the "hand shaking" is eliminate-d . It’s 
drawback is that the output is not truely real time because 
of the asyncronous operation of the AGT-10's. However, the 
reduced sampled data time and computer iteration time would 
make it a closer approximation of real time than presently 
exists . 
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APPENDIX E 

RTS5D MODIFIED PROGRAM NOMENCLATURE 



All 


A/D trunk 500 voltage 




A21 


A/D trunk 501 voltage 




A22 


Added mass coefficient in yaw moment equation 


A31 


A/D trunk 502 voltage 




A33 


Added mass coefficient in roll moment 


equation 


A34 


Added mass coefficient in pitch moment 


equation 


A4l_ 


A/D trunk 503 voltage 




ABSYD 


Y 

0 




ABSYOD 


V 




AD 


array of A/D lines 




A1COA 


speed line index counter 




AM 


mass 




APITCH 


pitch angle in degrees 




APRINT 


; f iterations between print execute commands 


ARR 


Real IARR (l) 




ARATE 


turn rate in degrees 




A ROLL 


roll angle 




AST OP 


performance index limit 




AXST 


initial condition X-coordinate on sea 


track 


AYST 


initial condition Y-coordinate on sea 


track 


BETA 


drift angle 




BOAT 


subroutine to generate perspective of 


SES 
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CASST 


automatic speed control mode indicator 


GDP ' 


bow real forces lumped coefficient 


CDX (i) 


drag forces lumped coefficients 


GDY 


sway forces lumped coefficient 


GDZP 


sidewall rolling moment lumped ceofficient 


DAL 


digital to analog call 


DELT 


iteration loop time 


DELTA 


commanded iteration loop time 


DFCCM 


total change in fuel command from speed comtroller 


DEC CM 1 


change in fuel command due to proportional speed 



controller 



DFCCM 2 


change in fuel command due to proportional speed 



controller 



DTI ME 


v time 


DTIMPLT 


scaled value for v 


DUCCM 


difference between ucom and u 


DXPLCT 


same as DTIMPLT 


DYPLCT 


scaled value for v 


DYREPET 


restart value for v 


SFLG 


error flag for TEXTC call 


ETIME 


r time 


EYREPET 


restart value for r 


SXPLOT 


scaled r time 


EYPLCT 


scaled r 


FI 


stem buoyancy force 


F2 


bow seal pitch force 


?3 


bow buoyance force 
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?C CM 


fuel command 


PS? 


port sidewall buoyancy force 


F5S 


starboard sidewall buoyancy force 


FTIME 


v time 


FTIMPLT 


scaled v time 


FXPLCT 


same as FTIMPLT 


FYPLCT 


• 

scaled v 


PYRE PET 


• 

restart value for v 


G 


gravitational acceleration 


GRAPKO 

GTIME 


subroutine to project image (non alpha-numeric) 
R time 


GTIMPLT 


• 

scaled R time 


GXFLOT 


same as GTIMPLT 


GYPICT 


• 

scaled R 


GYREPET 


• 

restart R value 


H (i ) 


dynamics maximum and minimum array 


HEAD 


craft heading 


HOLD 

HTIME 


subroutine to freeze display/program 

• 

u time 


HTIMPLT 


• 

scaled u time 


HYREPET 


analog restart y coordinate on sea track 


HYSESX 


analog y position of craft 


HYSSSY 


analog x position of craft 


IARR(l) 


digitized mode selector 


IARR(2) 


digitized effector angle 


IARR ( 3 ) 


digitized T7 
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IARR (4 ) 
IARR ( 5 ) 
IARR (6) 
IARR (? ) 
ICO 

I CCA 

I CRAFT 

ISR 

II?R 

I I PRD 
IIPRE 
IIPRF 
IIPRG 
IIPRH 
IIPRC 
IIPRP 
IIPRQ 
IIPRR 
IIPRS 
I IX 
IJ 
IJA 
IJAD 
IJAS 
I JAF 
I JAG 
IJAH 



digitized T8 
digitized T9 
digitized TlO 

digitized surge velocity command 
multiplex graphics index 
speed line index 
vehicle designator 
error flag 



time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 


time 


history 


limit 



print counter 
craft perspective counter 
craft perspective counter 
time history counter 
time history counter 
time history counter 
time history counter 
time history counter 
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IJAO 


time 


history 


counter 


IJAP 


time 


history 


counter 


IJAQ 


time 


history 


counter 


I JAR 


time 


history 


counter 


IJAS 


time 


history 


counter 


IJD 


time 


history 


counter 


IJE 


time 


history 


counter 


IJF 


time 


history 


counter 


IJG 


time 


history 


counter 


IJH 


time 


history 


counter 


IJO 


time 


history 


counter 


IJP 


time 


history 


counter 


IJQ 


time 


history 


counter 


IJR 


time 


history 


counter 


IJS 


time 


history 


counter 



IL 


graphics 


console 


display 


array 


ILAl 


graphics 


console 


display 


array 


ILA2 


graphics 


console 


display 


array 


ILA3 


graphics 


console 


display 


array 


ILA4 


graphics 


console 


display 


array 


ILA5 


graphics 


console 


display 


array 


ILA6 


graphics 


console 


display 


array 


ILA7 


graphics 


console 


display 


array 


ILA8 


graphics 


console 


display 


array 


ILA9 


graphics 


console 


display 


array 


ILA10 


graphics 


console 


display 


array 
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ILA11 


graphics 


console 


display 


array 


I LAI 2 


graphics 


console 


display 


array 


I LAI 3 


graphics 


console 


display 


array 


ILA14 


graphics 


console 


display 


array 


I LAIS 


graphics 


console 


display 


array 


I LAIS 


graphics 


console 


display 


array 


I LAI 7 


graphics 


console 


display 


array 


I LAI 8 


graphics 


console 


display 


array 


I LA 19 


graphics 


console 


display 


array 


ILA2 0 


graphics 


console 


display 


array 


ILA21 


graphics 


console 


display 


array 


ILAZ2 


graphics 


console 


display 


array 


ILA23 


graphics 


console 


display 


array 


ILA24 


graphics 


console 


dispLay 


array 


ILA2S 


graphics 


console 


display 


array 


ILA2S 


graphi cs 


console 


display 


array 


ILA2 7 


graphic s 


con sole 


display 


array 


ILA2 8 


graphics 


con sole 


displ ay 


array 


ILA29 


grapftics 


console 


display 


array 


ILAP 


graphics 


console 


displ ay 


array 


IPACK 


subroutine to load display arrays 


I PLOT 


graphi es 


console 


display 


array 



IPR time history counter limit 

IPRD time history counter limit 
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IPRE 


time history counter 


limit 


IPRF 


time history counter 


limit 


IPRG 


time history counter 


limit 


IPRH 


time history counter 


limit 


IPRO 


time history counter 


limit 


IPRP 


time history counter 


limit 


IPRQ 


time history counter 


limit 


IPRR 


time history counter 


limit 


IPRS 


time history counter 


limit 


I SWA 


display array 




ITEXT 


subroutine to display alpha-numeric data 


IVEIWA 


craft perspective/north road array 


IVIEWAA 


sea track array 




IVIEWB 


sea track axis array 




IVIEWC 


u time history array 




IVIEWD 


v time history array 




IVIEWE 


r time history array 




IVIEWF 


v time history array 




IVIEWG 


• 

r time history array 




IVIEWH 


u time history array 




IVIEWO 


p time history array 




IVIEWP 


q time history array 




IVIEWQ 


p time history array 




IVIEWR 


q time history array 




IVIEWS 


analog sea track array 


IVIEWZ 


time history axis array 
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IX 


x axis moment of inertia 


J 


time history counter 


JD 


time history counter 


JE 


time history counter 


JF 


time histroy counter 


JG 


time history counter 


JH 


time history counter 


JO 


time history counter 


J? 


time history counter 


JQ 


time history counter 


JR 


time history counter 


JS 

x. 


time history counter 
proportional speed controller gain 


KK 


intregal speed controller gain 


L3 


lever arm for buoyancy force 


LD 


draft of sidewall 


LP 


side thrust lever arm 


LZ 

N 


lever arm for sway drag force 
clock interrupt count 


NAD 


AD trunk line array 


ONEMOT 


subroutine to calculate SES dynamics 


CO 


lever arm for thrust side component 


oTir/is 


p time 


OTir.IPLT 


• 

scaled p time 


CXPLOT 


same as OTIIJPLT 


OYPLCT 


• 

scaled p 


CYREPET 


» 

restart p value 
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D 


roll rate 


P2 


plenum pressure 


FDCT 


roll acceleration 


PDOTM 


max roll acceleration 


PHI 


roll angle 


PMAX 


max roll rate 


PSI 


craft heading 


PTI-ffi 


• 

q time 


PTIMPLT 


• 

scaled q time 


PXPLCT 


same as PTIIYPLT 


PYPLCT 


• 

scaled q 


PYRSFET 


• 

restart q value 


Q 


pitch rate 


QDCT 


pitch acceleration 


QDCTM 


• 

max q 


QMAX 


max q 


QTir.E 


p time 


QTIMFLT 


scaled p time 


QXPLCT 


same as QTIMPLT 


QYPLCT 


scaled p value 


QYRSPST 


restart p value 


R 


craft turn rate r 


RDOT 


• 

r 


RDCTI/I 


♦ 

max r 


READCLOCK 


subroutine to sample real time clock 
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RESET 


subroutine to reset analog computer 


RHO 


density of water 


RMAX 


maximum yaw rate of craft 


RTIIS 


q rate 


RTIKPLT 


scaled q time 


RX 


I x 


RXPLCT 


same as RTII.TPLT 


RY 


J Y 


RYPLOT 


scaled q value 


RYRSPET 


restart q value 


RZ 




SSAST 


sea state 


SESX 


scaled X„ 
0 


SESY 


scaled Y 

o 


SF1 


scale factor in time history arrays 


SF14 


scale factor in time history arrays 


SLIP 


Beta influence coefficient 


S PD LIN 


speedline lower screen limit 


SPEED 


Craft's velocity in knots 



STARTCLCCK subroutine to start clock 



STIJ/TE 


analog time 


STIT/IPLT 


scaled analog time 


SXPLCT 


KYSESX 


SYPLCT 


HYSSSY 


SXREPET 


restart value for HYSSSX 


SYREPET 


restart value for HYSSSY 
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m 


total thrust 


T10 


thruster number 4 


TIOMAX 


max T10 value 


T7 


thruster number 1 


T7T-TAS 


max T? value 


T8 


thruster number 2 


T8KAX 


max value of T8 


T9 


thruster number 3 


T9 MAX 


max value of T9 


TEXTO 


subroutine to display alpha-numeric data 


TDCT 


T 


THETA 

TIM 


pitch angle 
u time 


Tins 


time 


TIMPLT 


scaled TIM 


TINT 


real clock interrupt frequency 


TITLE 


program name display array 


TITIEO 


program name display array 


TITLE1 


program name display array 


TM 


manual thrust command 


TMAX 


max T 


TSIDE 


sum of side thrust components 


tyaw 


sum of yaw moment 


u 


forward velocity (surge) 


UCOM 


forv/ard velocity (surge) command 
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UDOT 


u 


UDOTM 


max u 


um 


max u 


U-I.IUL 


scale factor in speedline generation 


V 


side velocity (sway) 


VCD 

VDCT 


thruster console pot array 

• 

V 


VDOTM 


• 

max v 


VMAX 


max v 


VS 


total velocity 



heave velocity 



WDCT 


heave acceleration 


VIE 


width of bow seal 



V/RITECLCCK subroutine to assign clock interrupts to 





variable 


WVJ 


lever arm for sway drag forces 


X 


Y o 


XC 


X 

0 


XCDOT 


Y 

' k o 


XDX 


boat perspective scale factor 


XGX 


boat perspective scale factor 


XPLCT 


display array 


X RE PET 


reset value for SESX 


XST 


scale factor for analog plot 


YO 


V 

"o 


YCDCT 


♦ 

V 

0 
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v*JT YY 


point D vertical position c on s tan 


q.t 


disolay array 




reset value for SZ3Y 


7QT' 

_ w — 


scale factor for analog plot 


z 


operational effector angle 


Z2 


upper effector angle limitation 


23 


lower effector angle limitation 


Zk 


zero angle upper limit 


25 


zero an ale lower limit 


ZA3 


scale factor for ZlO 


ZBRC 


effector breach limit angle 


ZI 


modified Z 


ZZZZ 


fixed step effector angle input 
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Appendix D 

RTS5D Modified Wiring Diagram 
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